Preparation of decoherence-free, subradiant states in a cavity 
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The cause of decoherence in a quantum system can be traced back to tlie interaction with the 
environment. As it has been pointed out first by Dicke, in a system of A*' two-level atoms where 
each of the atoms is individually dipole coupled to the environment, there are collective, subradiant 
states, that have no dipole coupling to photon modes, and therefore they are expected to decay 
slower. This property also implies that these type of states, which form an A — 1 dimensional 
subspace of the atomic subsytem, also decohere slower. We propose a scheme which will create such 
states. First the two-level atoms are placed in a strongly detuned cavity and one of the atoms, 
called the control atom is excited. The time evolution of the coupled atom-cavity system leads to an 
T— I ' appropriately entangled state of the atoms. By applying subsequent laser pulses at a well defined 

I time instant, it is possible to drive the atomic state into the subradiant, i. e., decoherence free 

■ subspace. Up to a certain average number of the photons, the result is independent of the state 

(N . 

of the cavity. The analysis of the conditions shows that this scheme is feasible with present day 
^ , techniques achieved in atom cavity interaction experiments. 

^ ■ 

1—5 ■ PACS: 42. 50. Fx, 03.67.Lx 

, Subradiant states of a system of two-level atoms has recently gained wide attention because of their excep- 

' tionally slow decoherence ■ This stability of quantum superpositions inside the subradiant subspaces originates 

^ , from the low probability of photon emission, which means very weak interaction between the atoms and their en- 

' vironment. Hence the subradiant states span a decoherence-free subspace (DPS) |p[-pl| of the atomic Hilbert-space 

' and consequently can become important from the viewpoint of quantum computation (QC) |10 1^. The scheme we 

o 



propose can be used to prepare subradiant states in a cavity. Our method is based on second order perturbation 
theory but the exact results verify the validity of the perturbative approach. We also investigate to what extent our 
scheme is independent of the state of the cavity field. Finally the requirements needed to prepare subradiant states 
in the proposed way will be compared with available experimental techniques ]T3|-|T6t . 

We investigate a system of N identical two-level atoms in a single mode cavity. Each individual atom is equivalent 
^ • to a spin- 1/2 system, and the whole atomic ensemble can be described by the aid of collective atomic operators J+, 
I J_ and Jz obeying the same algebra as the usual angular momentum operators . We consider the following model 
^ Hamiltonian: 

^ , H = Hf) + Hint = TiWa Jz + hojca'a + fig ia' + aJ+) , (1) 

^ ' where a and are the annihilation and creation operators of the cavity mode, Wa is the transition frequency between 
• ^ ■ the two atomic energy levels, lOc denotes the frequency of the cavity mode, different from Wa, and g is the coupling 
constant. We note that the Hamiltonian is written in the framework of Dicke's theory, i.e., with the assumption 
that all the atoms are subjected to the same field, which is a good approximation when the size of the atomic sample 
is small compared to the wavelength of the cavity mode. As discussed later in detail, there are experimental situations 
where this requirement is fulfilled. Our proposed scheme for preparing subradiant states involves a detuned cavity. 
We shall assume that the detuning is much larger than the resonant Rabi frequency: 

- = A > g. (2) 

Now any state of the atomic system and the cavity field can be expanded as a linear combination of eigenstates of 
i?o- These are tensorial products of collective atomic states and number states of the field: |j, m, A) ® |n), where the 
indices j, m and A label the atomic state (also called Dicke state [Q) while n refers to the nth Pock state of the mode. 
The quantum number j corresponds to the eigenvalues of the operator = + {J+J- -t- J_ J+)/2. This index is 



in one-to-one correspondence with the Young digram |17| that describes the permutation symmetry of the state. The 
possible values of j is A^/2, A^/2 — 1, . . ., the smallest value being if is even and 1/2 if A^ is odd. The index m 
of the |j, TO, A) Dicke state labels the eigenstates of the collective atomic operator J^, that is essentially proportional 
to the energy of the atomic subsytem. This is the index that is decreased (increased) by one under the action of the 
operator J_ (J+): 

J_|j, m. A) = y/ j{i + 1) - m(TO - l)|j, TO - 1, A), (3) 



1 



including the case when m = —j, when the result is the zero vector. The states with m = —j are the lowest ones 
of the Dicke ladders 0, they are called subradiant, because they have no diploe coupling to other lower lying states. 
Finally the index A distinguishes states with the same j and to. For more details see Refs. 0J|,|l^ 

Besides the collective atomic states \j,m,X), we shall also use the natural basis that assigns a well defined state 
to each individual atom. These vectors will be labeled by a string of 0-s and 1-s corresponding to the ground and 

12 N 

excited sates, respectively. E. g., the ground state of the atomic subsystem is written in this basis as | 00 ... O); this 
state (as well as the fully excited one) is also an element of the Dicke basis, |00 . . . 0) = \j = N/2, to ~ —N/2, A = 1). 

The form of the Hamiltonian (|l|) implies that the time evolution of the system shall exhibit two time scales: The 
first characteristic time is due to the self-Hamiltonian Hq and is approximately 27r/wa (or 2tt/uJc) and the second is 
proportional to 2t: / g. Generally g <C Wa ~ and the faster process induced by Hq can be eliminated by going into 
an interaction picture. However, if the frequency difference A is large enough, then the energy transfer between two 
adjacent eigenstates of Hq, differing in only one photon number, becomes negligible. This means that the amplitude 
of the corresponding collective Rabi oscillations will be very small, that is, the process on the second time scale will 
be unnoticeable and even slower mechanisms will become apparent. The situation is similar to the proposals and 

il- 

Hereafter we shall focus on the solution of the Schrodinger equation in the case when just a single atom is excited 
at t = 0. This initial state can be prepared by starting from the state |00 . . . 0), and exciting one well defined control 
atom. This excitation can be achieved via a third much higher lying level, so that the wavelength of the addressing 
pulse allows to focus it on the desired target atom |Q . For the sake of simplicity we always consider the control atom 
as being the first, hence the initial state will be written as 

|(/,(0)) = |100...0)(8)|n-l). (4) 

In order to find the complete analytical solution of the Schrodinger equation induced by the Hamiltonian (|^), in 
principle one should calculate all the eigenvalues and the corresponding eigenstates of H . Athough this problem 
can be solved analytically psj , more insight is given by a simple perturbative approach. The exact nonperturbative 
numerical solution of the Schrodinger equation verifies that results obtained via perturbation theory yield excellent 
approximations. 

The state 

|l)^(^-^f;|0...0l0...0)^®|n-l) = 

= \j ^ iV/2,TO = ~N/2 + l,\^ 1) ® \n~l), (5) 

which is in the completely symmetric subspace, and the subradiant states: 

\i) = \ j 7V/2 - 1, TO = -N/2 + l,\ = i-l)®\n-l) (6) 

with i = 2, 3. ..TV, have the same unperturbed energy, they span the iV-fold degenerate eigensubspace of iJg corre- 
sponding to the eigenvalue E^{n) — h{nuJc — Nu!a/2) — TiA. 

It can be seen that first order degenerate perturbation theory is not giving any correction to the energy, because 
all the matrix elements of Hint between the states above vanish, the action of Hint on vectors \ j, to. A) (g) |n — 1) gives 
a linear combination of | j, to — 1, A) ® \n) and |j, m + 1, A) ® |n — 2) that are orthogonal to the states (||) and (^. In 
order to obtain nonzero energy corrections we have to make a second order degenerate perturbation calculation , 
and find the eigenvalues of the matrix: 

E {i\Hint\Tn){m\Hinf\k) 

where the sum runs over all eigenstates of Hq with eigenvalc E^!^ ^ E'^ (n) . The only nonvanishing energy corrections 
in second order are the following: 

SEi = h^{Nn-2N -2n + 2), 

5E, = 5Ei + hN^, i^2,3...N. (8) 

At this point we can formulate the requirements that assure the validity of the perturbation theory: the magnitude of 
SEi and 6Ei must be much smaller than ?i|A|, the minimum of the difference between E^{n) and all other unperturbed 
energy levels. 
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The most important consequence of Eqs. (g) is that the Bohr frequencies that determine the time dependences of 
the subradiant and non-subradiant states are different. 

Now we expand the initial state as the hnear combination of the fully symmetric (non-subradiant) state 
and an appropriate subradiant state: 



|2> 



1 



(AT - 1)|100...0) 



N 



^|O...OlO...O) 



(9) 



By assigning the symbol |2) to the state in Eq. 
subspace. Now the initial state reads 



we have utilized the freedom of choosing a basis in the subradiant 



I<^(0)) = -^|1) 



N -1 
N 



|2>. 



(10) 



By the aid of this expansion and using the Bohr frequencies resulting from (H), it is easy to calculate the time evolution 
of the state (p^. Discarding an overall phase factor, this time dependent state has the form 



\m 



^^exp (iN^t]\l) + 



N ~1 
N 



|2), 



(11) 



or, on using Eqs. (^ and (||): 



|(/,(t)) = {N cos{at) - i{N - 2) sm{at)) \100 .. .0) 
+2i sin(at) ^ |0 ... 010 ... 0) (g) |n - l)/iV. 



k=2 



Here we introduced the parameter 



2A ' 



(12) 



(13) 



which is independent of n. Because of this latter fact, from now on the state of the cavity field will be omitted in the 
notation. We also note that the characteristic time of the time evolution, 2TT/a, is much longer than that of the free 
evolution due to Hq, being the consequence of the fact that the evolution described in Eq. ( [ll| ) is induced by a weak, 
nonresonant interaction. 

Eq. ( |l^ reveals that in \4>{t)) the weight of the state |100 ... 0) and those of the states with the first atom unexcited 
changes during the course of time. As we can see, the moduli of the corresponding coefficients in Eq. (O) are 



iV2 cos2(at) + {N- 2)2 sin2(at) 21 sm( at)\ 
and 



N 



N 



respectively. Comparing these values to Eq (^, it can be shown that for an arbitrary N there exists a time instant 
tm when 



|(/.(t)) = , ^ \{N - \)e'^ 1 100 ... 0) 



N 



^|0...010...0) 



k=2 



(14) 



which differs from the subradiant state |2) only in the phase factor e*"^ of the first term. Combination of the previous 
two equations and Eq. (^ yields the following requirement for tm- 



N"^ cos^ {atm) + iN-2y sm^iat,^) 



|2sin(ai„)| 



= N-1. 



(15) 
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We can find a solution of this equation for all > 1: 



Sinai™ = ^N/{AN-A), (16) 

and also obtain cos (p — 21^-2 (0 ■ 

Now it is clear that at the time instant given by Eq (^6|), an appropriate rapid change in the phase of the state 

k 

1 100 ... 0) relatively to all other states |0 ... 1 ... 0) leads to the subradiant state |2). 

On the other hand, Eq. ( p^ ) also shows that the required phase transformation is equivalent to the elimination of 
the phase difference between the |l)c excited and |0)c ground state of the control atom. Therefore we consider the 
action of a strong laser pulse on the control atom. In order to obtain precise addressing ||2^], the laser is to be tuned 
in resonance with an allowed transition |l)c — > |e)c, where |e)c denotes a state of the control atom with much higher 
energy than \l)c- E. g., by the appropriate choice of the phase of the complex Rabi frequencies of two tt pulses leads 
to the phase transformation reqired to prepare the subradiant state |2). Additionally, the duration of a Rabi period 
due to the strong, resonant laser pulse is much shorter than the characteristic time that governs the time evolution 
(^). We note that the idea of introducing phase transformation in a multilevel system by the aid of short laser pulses 
has appeared in a somewhat different context in p7[ |. 

Now we show that our scheme is independent of the state of the cavity field, and write more generally the initial 
state as: 

|(/)(0)) = I100...0)® IV'W) = I100...0) ® ^c„(i)|n). (17) 

n 

We use the fact that the interaction Hamiltonian Hint does not mix states with different number of excitation 
(essentially n + m): 

{j,m,X\® {n\H,nt\n') ®\j,m±l,\) =0, (18) 

unless n' = n =p 1. This implies that the calculations based on second order perturbation theory can be performed 
for each A^-fold degenerate energy level of iJo corresponding to different values of n. After replacing the state \n — 1) 
with IV'(i)) in Eqs. (§ and (|), we obtain the following result: 

5E,^5Ei = ^^NY,\cn\'' ^2a, (19) 

n 

which is therefore also valid in this general case. Thus we have proven that our scheme does not require special 
preparation of the cavity field. However, it should be borne in mind, that the results above are based on perturbation 
theory. For given TV, g and A the validity of the perturbative calculations depends on (n) , the average number of 
photons in the cavity field. Hence it is clear tha t our scheme can not be independent of the average photon number 
on a very large scale. Nevertheless, until ■^yjN{n) <C 1, all the previous statements hold. For the case of 10 
atoms, we have performed exact (nonperturbative) numerical calculations for the experimentally realizable ratio |lj] 
oi g/ A = 30, and found that the time evolution follows Eq. (^l|) within 2% relative error in the coefficients. 

Finally we compare the requirements of our scheme with the experimental possibilities. The atom cavity experiments 
of Haroche and co-workers, as described in the review paper [ p8[ , show that the description of the interaction of a 
number of Rydberg atoms with a single mode cavity is truly described in the framework of the Dicke model. In 
more recent experiments [p^|l6| the interaction of a detuned cavity with one and two atoms has been found to be in 
agreement with theoretical predictions. The parameters realized in these experiments with rubidium Rydberg atoms: 
17/277 « 2AkHz and detunings as large as A/27r w SOOfciJz show that the conditions of the validity of our perturbation 
approach hold as much as for about hundred atoms, because the average photon number in a cavity can be kept much 
less than 1. Taking for instance 10 atoms and g/ IS. — 30, we have a = 2.5 x 10'*l/s giving for tm a value of 22/is. The 
interaction time of the atoms and the cavity must be longer than t,-a what can be achieved already with atoms with 
somewhat less than thermal velocities for centimeter sized cavities. Finally we note that addressing of single atoms is 
a common problem in almost all of the proposals in QC, but there are promising works indicating future success [p4| . 

In conclusion, we have proposed a method to prepare decoherence-free, subradiant states of a multiatomic system. 
We have shown that our perturbative approach is compatible with present day techniques in atom cavity experiments. 

We thank S. B. Zheng for discussions. This work was supported by the Hungarian Scientific Research Fund (OTKA) 
under contract T32920, and by the Hungarian Ministry of Education under contract FKFP 099/2001. 
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